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FIELD OF THE INVENTION 



The invention relates to tape systems, and in particular, to systems for 
determining the tape pack radius of supply and/or take-up reels. 



BACKGROUND OF THE INVENTION 

Tape systems generally have two reels for storing tape, namely, a supply reel 
and a take-up reel, a capstan for moving the tape from reel to reel and tension arms for 
regulating the tape tension. High-performance tape systems also include servo 
systems, which regulate tape position and velocity. The servo systems rely on 
+10 estimates of the tape pack radii to determine how to control the rotational speeds of the 
reels to achieve the desired tape velocity and position. The more accurate the 
estimates, the more precisely the servo system can control the movement of the tape. 



Ui 



all aspects of the system operations. For example, good estimates are important in 
15 determining from which reel to draw the tape to wrap around a scanne^An inaccurate 
estimate could result in an over-rotation of the selected reel. ^Farther, good estimates 
are important to determining when to decelerate a higrT-speed rewind operation, again 
to avoid over-rotation of one of the reels thafmay result in the breaking of the tape of 
the detatchment of the tape fronrpHe reel. Also, good estimates are important to 
20 determine if there is sufficient tape available on the supply reel to complete a record 
c p e ra t i c rv^-ffta ecu rate estimates of the reel pack radii can result in incomplete record 
op^yeftsriniie syslUm sulliclUIIUy unuei es lilies II IU type 1 pmroa. 
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In prior known systems the tape pack radius is calculated from measurements of 
the angular positions of the reels and the capstan. The position measurements are 
made by, for example, optical encoders that count the number of slots that pass 
between a photo detector and a light source as the reel rotates. The calculations 
5 produce results that are at best as accurate as the position measurements, which tend 
to be "noisy." With optical encoders, for example, the measurement noise is due in 
large part to quantitization errors. At slow speeds these systems tend to produce 
relatively inaccurate results because the position measurements are comparable to the 

i 

f 

! quantitization errors. 



filters, in order to smooth them. However, this approach has two significant problems. 
First, the signals produced by these filters always lag behind the true tape pack radii, or 
in other words, the estimates are biased. Second, these filters are slow to converge. 
Moreover, there is an intrinsic tradeoff - the more the filter smooths the output signals, 



15 i.e., the estimates of the tape pack radii, the more lag is introduced into the system and 
the slower the convergence. 

SUMMARY OF INVENTION 

The invention is a tape system that uses a Kalman filter to produce both an 
estimate of the tape pack radius and an estimate of the associated estimation error 
20 variance. The Kalman filter produces the estimates based on current and prior angular 



estimates of both the tape pack radius and the estimation error variance. The filter 
gain is based on the prior estimation error variance and a calculated measurement 



ylo 



Certain prior systems have processed the noisy measurement using low-pass 




isssuroments of the reel, capstan, and tension arm. and the prior 
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error variance, which is a function of the current and prior measurements. Preferably, 
the system uses two Kalman filters, one to produce the estimates for the supply reel 
and one or produce the estimates for the take-up reel. 

More specifically, the Kalman filter includes a model of the dynamics of the 
5 system that it uses to predict what the tape pack radius and estimation error variance 
will be at the time of a new measurement. The Kalman filter then updates the estimate 
of the tape pack radius when a new measurement is made, based on this predicted 
radius and a correction factor that is the weighted difference between a "measured" 
reel pack radius (calculated from reel, capstan, and tension arm measurements) and 
Jo the predicted radius. This weight, which is Kalman filter gain, is a function of the 



tfB BRIEF DESCRIPTION OF THE DRAWINGS 

The above and further advantages of the invention may be better understood by 
referring to the following description in conjunction with the accompanying drawings, in 
which: 



Fig. 3 is a flowchart of the operations of the system of Fig. 1 when the system is 
operating in a coasting mode. 



estimation error variance and a measurement noise, or error variance. 



Fig. 1 is a functional block diagram of a system constructed in accordance with 



20 the invention; 
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DETAILED DESCRIPTION 

1 . The System 

Referring to Fig. 1, a tape system includes a capstan 14 that transfers tape 12 
between a supply reel 10 and a take-up reel 16. When the reels rotate clockwise, as 
5 depicted by the arrows, tape is wound off of the supply reel and onto the take-up reel 
and, thus, the tape pack radius of the supply reel decreases as the tape pack radius of 
the take-up reel increases. 

Angular position transducers 21-25, which may be, for example, optica! 
* encoders, are attached to the shafts (not shown) of the capstan 14, the two reels 1 0 
1=10 and 16, and the tension arms 1 8 and 20. The transducers operate in a conventional 
manner to measure, respectively, the angular positions of the capstan, the reels and 



A supply Kalman filter 26a receives measurements from transducers 21 , 22 and 
24, and a take-up Kalman filter 26b receives measurements from transducers 23, 22, 



vfll5 and 25. These measurements are used by the filters to produce estimates of the tape 
pack radii and estimates of associated estimation error variances, which provide a 
confidence interval around each tape pack radius estimate. The operations of the 
filters 26a and 26b are discussed in more detail with reference to Fig. 2 below. Also as 
discussed in more detail below, the estimates produced by the filters are used by a 
20 servo system 28 that controls the motion of and the tension in the tape 12, and by a 
tape length estimator 30 that determines how much tape is available for a record 
operation. 



gfaHdes r ^e anners T Jonaitudinal heads. and-SO forth, no ne of w h i fi h.a co ohow n io Jhe 



the tension arms. 
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dfu w i n g. AfrOowg qo those o l oroo a to oro - tixad I h position MU that they do not causelfe 
tape path to vary, their presence is immaterial to the operations ofth^JSaJmaTTfilters 
26a and 26b. The effects on the operaiipjas-eHheTnters^ a change in the tape path 
lengJMorexSmple, by means of movable guides, is discussed below with reference to 

2. The Theory 

We can look at the system as two subsystems, namely, a supply subsystem that 
includes the supply reel 10, the supply tension arm 18 and the^capstan 14, and a take- 
up subsystem that includes the take-up reel 16, the take-up tension arm 20 and the 
capstan 14. The two subsystems have identical theories of operation, and we discuss 
herein the theory of operation of the take-up subsystem. 

Let us consider the case when position measurements are made at regular short 
intervals of, for example, 20 milliseconds. If the tape is moving, the length of tape that 
is wound onto the take-up reel 16 must travel past the capstan 14, and we can state 
that, to a high degree of accuracy, 

r^0 r + MAO a -r c A$ e = O 
where r r is the radius of the take-up reel tape pack, A9 r is the change in angular 
position of the take-up reel 16 from the previous sample time to the current sample 
time, A9 C is the change in angular position of the capstan, u a is a gain determined by 
the geometry of the tape path and the nominal position of the take-up tension arm 20, 
and A6 a is the change in angular position of the arm. The radius of the take-up reel is 
thus: 
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r c A0 c - n a L0 a A 
r ' = \a eC ^ n - 1 



If the tape is moving relatively slowly the "deltas," that is, the changes in the 
angular positions are comparable in magnitude to the position measurement errors. If 
the system relies solely on these measurements to determine tape pack radius and 
5 ignores both past measurements and the dynamics of the system, errors in position 
measurements translate directly to inaccuracies in the determination of the tape pack 
radius. The prior systems discussed above in the Background section all suffer from 
Cjj this problem. 

% Our system circumvents those problems through the use of a Kalman filter, 

JjiO which has as an input signal a "measured" radius r,„, that is calculated using equation 

yj 

R 10. The Kalman filter incorporates a model of how the tape pack radius changes as a 

m 

]_ u function of reel position to predict an estimated tape pack radius. The filter then uses 

£ the measured radius to correct this prediction. The derivation of the predictive model is 

m 

considered here. 

W 1'5 Consider the tape pack as a hollow, circular cylindrical mass with width V, 

and a uniform density "d," i.e. mass per unit volume, where d is the density of the tape 
itself. If a length of tape " / " units long, with width w and thickness "8 " is wound onto 
the reel, the mass of the tape pack increases by the mass of that length of tape, and 
the radius of the tape pack increases accordingly. 
20 The mass of the tape pack before the length of tape is wound on is: 

dwnr* '- dwnr 0 2 

where r f is the outer radius of the tape pack and r 0 is the inner radius of the tape pack, 
i.e., the radius of the reel hub. Suppose that the radius of the pack increases to r 2 
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when the length of tape is wound onto the tape pack. The mass of the tape pack thus 
increases to: 



The increase in the mass of the tape pack due to the tape being wound onto the reel is 
dwl<5> which is the mass of the additional tape, hence: 

dwnr 2 2 - dW7iro 2 = dwarf - dwnr 0 2 + dwIS. 

Solving for r 2 we have: 



Given an initial reel radius r, this equation enables us to approximately predict 
what the radius will be after we wind on the additional tape of length /. What we are 
really interested in, however, is an equation that allows us to predict what the radius 
will be given an initial reel radius and an amount by which the reel rotates. If the reel 
rotates in the clockwise direction by a small amount A0 r and the reel radius is initially r, 
a length of tape equal to rA0 r is wound onto the tape pack and the tape pack radius 
increases by Ar, which according to eqn. 3 satisfies: 



dwTor 2 - dw7ir 0 2 



eqn. 2 




eqn. 3 






hence 



8 



and, in the limit, 
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Ar = 



A0 r S 
2nr 




jveequation isthg 
the time tf for which the prediction is required: 



eqn. 4 



gral ofeqii. 4, lium ail initial UkiAHv to 



r(tf)-r(to) = ^[e r (t f )-e r (t 0 )^ 



In 

where r(tj) is the tape pack radius at time tj, apd^ r (tj) is the angular position of the take- 
up reel at time t|. Thus: 

Inaccuracies in this mocteffare accounted for in the Kalman filter by an estimate of 
"noise" applied toiHesystem - that is, by quantifying as applied noise essentially 
unmeasuraWeattributes of the system dynamics. The applied noise is one of the 
factors/tnat is used in determining the Kalman filter gain, k, as discussed below with 
reference id hiys: 2 H aM<4-3v-» 

To determine the filter gain, the system also requires an estimate of the 
measurement noise, or error, variance. The system implicitly assumes that the errors 
in the various position measurements can be represented by independent, zero-mean 
random variables 8j, which each have "small" variances. We know that for a smooth 
(differentiable) function g that maps V real numbers to one real number, we have: 

var(g(x 1 +8 1l ... l x n +sn))« 
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varCO 



From this relation we obtain the approximate measurement noise variance. Based on 
equation 1 , we define the function g as: 



Q(0c, Oe-,9; 0^., 6 r , Of.)- 



0 -0 



I where 6 C is the position of the capstan at the current sample time, 6c is the position of 

l the capstan at the previous sample time, 0 a and 6a- are the positions of the tension arm 

J at the current and previous sample times, and 0 r anc j Q r are the positions of the take-up 

fo reel at the current and previous sample times. If these quantities are perturbed, 

i- respectively, by independent, zero-mean random variables e c , e a and e r , the 

v approximate measurement error variance o m 2 is: 

T 

G » = ^ £ c) + M 0 2 ^rie a ) + g(O c> O ct 0 a ,O a .AA-) 2 vai(*,)] eqn. 6 



15 The variances var(e e ), var (e a ) and var (e r ) are computed by assuming that the 

quantitization error is a random variable with zero mean that is uniformly distributed 

2/r 

over an interval of size — radians, where N is the number of encoder counts per 

N 

revolution. 

gjy^ Tfa a^tertLifl Bert^ tii^bUlenibiilb-thrtHKJtJltfi mines to be uiiiedsuiidbfe 
20 discussed below with reference to Fig. 3, the systi 

unreasonabJ&jHhe associated measurement error variance, am, is greater than the 
maximum varia nce iw 7 , u i if the l l neg-b i y i na i iie asui e ii ie ii rn t iw val arouna tne 



rs a measui einei it to be 
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qateotelg^jaeasafeU r ad iu s i' m , dsTJelei inliie d b y the m e as u r e 
least partially included within the interval fromjTejnwtt^M tape pack radius, r^n, to the 
maximum tapejDacJsja^HasTFmax. Ignoring the unreasonable measurements greatly 

Ttfe Kalman filter to capstan snpr 



imprc 



tness 



The Kalman filters use equations 1 and 6 to produce estimates of the tape pack 
radii and the associated measurement error variances, respectively, as discussed 
below. 



4t0 



tfl5 



3. System Operation 



ot>ortaporfei 3 '0> 



20 



wound on the supply and take-up reels 10 and 16, is loaded into the system, 
system has no indication of the tape pack radii. Accordingly, the taeetJack radii 
processor 26 in step 300 must initialize the Kalman filter wjth^l) an initial tape pack 
radius estimate, (2) an initial estimation error varjaffce and (3) initial position 
measurements of the capstan, the tensjaffarm, and the reel, that is, e c .9 a ,and 6 r , where 
e r . is the angular position of tfje'fape reel under consideration, for example, the take-up 
reel. We discuss belpw^the operations of the system in determining the tape pack 
radius of theiaKe-up reel 16. The system performs the same operations to produce 
of the tape pack radius of the supply reel 10 and the associated estimation 



The Drocessor 26 uses as the initial estimate of the tape pack radius: 



r = 



'max mm 
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which is the radius when one-half of the tape is wound on the reel, and uses as an 
estimate of the initial estimation error variance: 



which is the variance of a uniformly distributed variable over the interval r min to w 
This initial estimate of the error variance is relatively large and indicates that the initial 
estimate of the tape pack radius is very likely to be inaccurate. 

The processor 26 takes a next set of position measurements G c .9 a . and G r , and 
calculates A9 Cl A9 a and A9 f , using the initial set of measurements as the measurements 
of the previous sample time, namely, 8c., 9a- and 8 r . (step 302). 

Using these delta values, the processor, in step 304, calculates a "measured" 
radius, r m : 



The processor, in step 306, next determines if the calculated radius falls outside of the 
interval of possible radii, namely, the interval r mi „ to r^. If so, the processor sets the 
calculated r m to r^n if r m is less than r mi „, or to r™ if r m is greater than w (step 308). 

The system then in step 310 calculates a measurement error variance o m 2 , using 
the measured A9 r and the calculated radius r m : 



max min 



12 



r c A0 c - Ma A0 a 
A0 r 



a « 2= T^-k 2var «>+ 
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If the calculated measurement error variance is larger than the maximum estimation 
error variance, which is determined by equation 7, the processor determines that the 
angular position measurements are unreasonable (steps 312, 314). 

Assuming the calculated measurement error variance is not greater than the 
maximum estimation error variance, the processor further tests the reasonableness of 
the measurements by determining if at least a portion of the three sigma interval around 
r m is within the interval of r^ to r™ (step 316). If so, the processor determines that the 
measurements are reasonable, and uses the calculated radius r m and measurement 
error variance o m 2 in the Kalman filter to determine estimates of the tape pack radius 
and the estimation error variance. Otherwise, the system determines that the 
measurements are unreasonable and, as discussed below, uses predicted values for 
the estimates. 
y^V^ M ^ h< ^iRi n l ?''r p m^n t T i are reasnoahl o. tho dk 
Kalman filter gain, k, as: 



where v* is the estimation error variance extrapolated from the previous sample time. 
The processor next, in step 320, updates the estimate of the tape pack radius as: 

and updates the estimation error variance as: 
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These updated estimates are the output signals of the filter. The servo 
controller 28 and a tape length counter (not shown) then use these estimates to 
determine how to control the tape system, as discussed below. Next, the tape pack 
radius and estimation error variance are propagated to the next sample time. The 
radius estimate is propagated according to the predictive equation: 

r 

which comes from equation 5, and the estimation error variance is propagated 
according to the equation: 
v~=y*+a v 2 , 

where ov 2 is a constant that represents the inaccuracies of, or the noise applied to, the 
predictive model. 

The servo controller 28 uses the estimated tape pack radii in conventional 
feedback and feed-forward loops that control the rotational speeds of the reels and 
move the tape at a desired velocity. The tape length estimator 30, which determines 
how much tape is available for record operations, uses as a lower bound for the tape 
pack radii: 

r-3a. 

This is a conservative estimate, which prevents the tape length estimator from over 
estimating the length of the tape available for a given record operation. 

If the processor 26 determines that the measurements are unreasonable, it does 
not update the filter gain. Further, it does not calculate the updated estimates. 
Instead, it predicts what the estimates should be based solely on the previous 
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D 



estimates. The processor, in step 315, sets f* = f~ and v* = v'and predicts an 
estimate of the tape pack radius r~ as: 

In r 

and the error variance v as : 

V"=V + +(T V 2 

The processor then sets the measured values 6 C , 8« and 9 r as the prior values e c ..6a-, 
and e r ., respectively, and repeats steps 302-320 to calculate updated estimates for the 
next sample time. 

The system continues to determine and, as appropriate, predict estimates for the 
jjo tape pack radii and the associated error variances as long as estimates are required, 
^ that is, as long as the reels continue to rotate. 

tfl tape path length is altered by anything but the tension arm. For example, the^stem 
* operates in the coasting mode when the tape is being unwound fron>dne reel and 
15 wrapped around the scanner (not shown). The system atep'dperates in a coasting 
mode when the tape is being withdrawn from the*s€anner. These operations; which 
significantly change the length of thejatfepath, are not included in the underlying 
assumptions upon which the^fteasurement model in the Kalman filter is based. 
Accordingly, the meaiured radius r m produced by that model cannot be used to 
20 determine tbe^estimates of the tape pack radius and the error variance when the path 
lengtK ^^ alte red^ he ^ stimato ^ ^ nrp fl ict i v o m nrlnl > 
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^Jj^Y T te , svst awHhas^xtEapolata^ o m t he current e sti ma te s o f t ap a p a a? rc 
estimation error variance to the next sample time by taking a next set of position 
measurements (step 402), determining A9 C , A9 a and A9 r (step 404) and estimating the 
tape pack radius by: 



r=r + +—A0 r 
In 



and the error variance by 



v =v +<r 



sift' 

■ si* 
HP 

w 



where ov 2 is a constant that represems the inaccuracies in the predictive model (step 
406). It next, in step408fsets 8 C ., 9a- and 9 r . equal to the measured values and f + = f 
and v + = v , aadat the next sample time again begins the process of predicting the 
estimates 



15 




20 



varying the tape patn lengtn. unce tnetaBfijjawrrengin is neld constant, the processor 
a g a ii%^f^ ^hfvftttnf -r|nin . nnd-o peiates^ disPi isspH ahnve^itlijeferenceJaBaJa. 

Generally, the system operates in the coasting mode for only a small number of 
consecutive sample times. By predicting new radii and estimation error variance 
values during the coasting operations, the processor can determine at the end of these 
operations reliable updated estimates for the tape pack radii and error variance, as well 
as a reliable filter gain, which is based in part on these estimates. 

raUlU b'e revised to include the" 




eceliSry'since the time the system speflfls 
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se-Httl a llt »e p erforming these operations. Thus, the trade oif orTncrerasetTcomplexji 
In the model for more accuracy during these limited pat h-varying opei attofts-seems 
unwarranted. Further, since the estimaJesifl^aYiced by the Kalman filter converge 
rapidly, the acajEaey-oTTfiesystem is only minimally reduced by not including in the 
moqei tno pnthvniTnTtr^P^rnTinrrr'''' 

The foregoing description has been limited to a specific embodiment of this 
invention. It will be apparent, however, that variations and modifications may be made 
to the invention, with the attainment of some or all of its advantages. Therefore, it is 
the object of the appended claims to cover all such variations and modifications as 
come within the true spirit and scope of the invention. 
What is claimed is: 
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